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Abstract 

A systematic expansion of the many-body master equation for active matter, in which motors 
power configurational changes as in the cytoskeleton, is shown to yield a description of the steady 
state and responses in terms of an effective temperature. The effective temperature depends on the 
susceptibihty of the motors and a Peclet number which measures their strength relative to thermal 
Brownian diffusion. The analytic prediction is shown to agree with previous numerical simulations 
and experiments. The mapping also establishes a description of aging in active matter that is also 
kinetically jammed. 
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At the nanoscale and above we encounter active forms of matter in which internal or 
external energy sources supplement passive Brownian motion to allow large-scale structural 
rearrangements. Constant agitation by biomolecular motors and force-generating polymer- 
izations allows the cytoskeleton of eukaryotes to undergo adaptive dynamical and structural 
changes in response to environmental disturbances. Other forms of active matter range from 
assemblies of entire microorganisms [l| to collections of artificial microscopic swimmers which 

nfi n 

are self-propelled |2|, |3[ and turbulently agitated active suspensions on the colloidal scale 

In this paper we show that, the steady-state statistical mechanics of these diverse out- 
of-equilibrium forms of active matter as well as their fluctuations and responses can be 
described using the concept of an effective temperature. The notion of effective temperature 
has been useful in describing passive glassy systems jsl, weakly driven systems such as gently 
sheared supercooled liquids and glasses js, 7| and vibrated granular matter j^, as well as in 
approximate theories and simulations of active biological matter jo- 12|. 

Modification of the fluctuation-dissipation theorem, signalling the nonequilibrium nature 
of the active processes like those in cells, has been directly observed in the mechanical prop- 
erties of reconstituted cytoskeletal networks [isl [l^. These measurements were consistent 
with an effective temperature nearly 100-fold higher than the ambient thermal temperature. 

In this Letter we will show that the effective temperature for systems with motor-driven 
active processes can be directly related to the properties of the motors. We describe the 
motors as generating a time series of isotropic kicks leading to a master equation description 
for the dynamics. By investigating this master equation in the small kick limit we obtain 
an expression for T^s explicitly depending on the total motor activity and the susceptibility 
of the motor dynamics to imposed mechanical forces. Motors lead to an enhanced diffusion 
for the active system regardless of susceptibility. Our theory predicts also a linear relation 
between the effective diffusion coefficient and the rescaled effective temperature, which de- 
pends on motor susceptibility. Our prediction is confirmed by recent numerical studies of 
systems with adamant motors 15|. 

An important consequence of our analysis is that the effective temperature allows a simple 
rescaling of the equilibrium phase diagram to be used for the nonequilibrium motorized 
system. We show that the effective temperature coincides with the results based on previous 
self-consistent-field calculations [igI. Active matter often is jammed and in a glassy state. 
The present analysis implies that Lubchenko-Wolynes aging theory 17| can be taken over 
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to active matter simply by assigning the ambient temperature in that theory to be the 
motorized effective temperature and treating the fictive temperature as a history-dependent 
dynamically controlled variable as in structural glasses below their glass transition. 

We start by following Shen and Wolynes [l^ who modeled the stochastic nature of motor 
kicking via a master equation for the many-body probability distribution function \l/({r},t) 

^vl/({f},t) = (LFP + LNE)*({rl,t). (1) 

Here Lpp = D^J^i^i ' ~ Dof3J2i^i ' (~Vjf/) is the usual many-body Fokker-Planck 
operator describing passive Brownian motion with Dq denoting the thermal diffusion con- 
stant at bath temperature T and (3 = l/ksT. The gradient of the many-body interac- 
tion free energy function U{{r}) = U{fi,r2, ■ ■ ■ ,fv) = X]<ij> ^(^« ~ "^j) gives the local 
force thus the thermal drift motion of individual particles and ri is the position of the zth 
particle. The effect of the nonequilibrium motor processes is summarized by an integral 
kernel L^^m{{r},t) = J Uid?i[K{{r'} ^ {f})^{{r'},t) - K{{f} {r'})$({f}, t)], where 
K{{r'} — !■ {r}) encodes the probability of transitions between different particle configura- 
tions. The (motor) kicking noise is thus a finite jump process with a rate that depends on 
whether the free energy is increased or decreased by a motor step 

k = K[e{AU) exp(-s„/3Af/) + 0(-Af/) exp{-Sd/3AU)]. (2) 

Here is the Heaviside step function and AU = U (^f+ fj — U (r) is the energy change due 
to the kick identified by a vector /. The kick step size / and the basal kicking rate k define 
the dimensionless motor activity A := kI'^/Dq which is analogous to the Peclet number in 
turbulent diffusion. This model rate couples the chemical reactions leading to the motor 
activity to the local mechanical forces acting on the motor; the assumed dependence upon 
instantaneous particle configuration reflects an assumed Markovian character of the dynam- 
ics without significant time delays. This is an idealization of the biochemical mechanism 
of real motors that doubtless possess intermediates in their function. We parametrize the 
coupling of a motor to the external forces it must overcome by the motor susceptibility s 
which may take different values for uphill moves and for downhill (s^) moves. s„ and Sd 
depend on the chemical mechanism of the motors. When s — )■ 1 the motors are susceptible; 
they slow down when they go up against mechanical obstacles and accelerate when they 
take energetically downhill steps; in contrast s — corresponds to completely adamant 
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motors which kick at an unperturbed rate being unaware of the instantaneous free energy 
landscape. Adamant motors use and waste a lot of energy. 
Explicitly we write 



LNE^({r},t) = /€^ jdnj dr[ | 5{n - r\ -I )w[U{- ■■,?,,■■■)- U{- ■■ ■■ ■)]^{{r'},t) 

- 5(r-;-r', + r)ti;[f/(--- ,rl,---)-f/(--- ,r'^,---p(H,^X 

The angular integration denoted by j dh averages over possible directions of kicking. For 
simplicity we provide results for the isotropic case. Again for simplicity we assume a fixed 
step length in the pair of delta functions. Our description of the rates gives w\Ui — U j] = 
e{Uf - U,)exp[-Suf3{Uf - U,)] + e(t/, - Uf)exp[-Sdf3{Uf - U^)]. 
When Su = Sd = s, one finds even more simply 



■sli[U{r,)-U{ri-l)\ 



-e-^^[^(^^'+'^)-^(^^^%({--- ,f,,---},t)|. (4) 



We assume the kicking of different particles at any time is uncorrelated. 

To obtain our promised results, we expand the distribution function and the kinetic rate in 
powers of / up to the quadratic order. This immediately leads to an effective Fokker-Planck 
equation 

-vl/({f},t) = Deff 5^ {V.^vl; - V. ■ [{-V.^sU)% }, (5) 



where 



^eff = I5o(l + i^!^), (6) 



2d a 



{MPV = = ( 1 + + ) . (7) 



2dDoJ / \ dDo 

Here the identity {cos^ 6)n = l/rf is used for general spatial dimension d. Rotational sym- 
metry in motor susceptibility eliminates odd powers in / upon the angular integration. 

These simple expressions ©-(ITj) valid for general dimensions have nontrivial implica- 
tions. First, in the small kick limit (kick step size is small compared to average particle 
separation yet the kicking rate can be quite high), the active system, while far from equi- 
librium in a strict sense, behaves as if it is at an effective equilibrium characterized by an 



effective temperature Tes which will however be quite different from the thermal tempera- 
ture T. It follows that the steady state is described by the effective Boltzmann distribution 
as \l/({r}) a exp[—U{{r})/kBTcs]. Second, the approach to effective equilibrium is governed 
by an effective diffusion constant i^cfr (Eq. E]) which is enhanced by the active processes re- 
gardless of the motor adamancy. This is consistent with recent experimental observation on 
enhanced cytoplasmic diffusion which plays such an important role in cytoskeletal assembly 
[l^ . Finally, it is clear from Eq. ([7]) that Tcs can be predicted knowing only the motor 
activity A = kI'^/Dq and the motor susceptibility s; susceptible motors with s > 1/2 yield 
Tcff < T; more interestingly, for very high motor activity, i.e., A ^ 1, the effective tempera- 
ture diverges as T^s/T ~ 1/ (2s) as s — > 0, indicating that intense kicking by adamant motors 
leads to a very high effective temperature just as observed in experiments and simulation 
studies. In principle motors could have slip bonds [2^ leading to negative s and negative 
effective temperatures, a situation we will explore in another publication. 

The results are easily generalized to the case of asymmetric susceptibility. It is not 
difficult to show that the factor (s^ — s^) accompanies all the cubic-and-above odd powers in 
/ and thus doesn't modify the effective equilibrium at quadratic order. Eqs. ©-(IT]) remain 
intact except for a direct substitution of {su + Sd)/2 for s. The effective temperature depends 
only on the sum of uphill and downhill susceptibility, and if s„ + = 1 (i.e. s = 1/2) we 
have Tcff = T. 

These predictions can be easily verified using results for various motile systems al- 
ready studied numerically or experimentally. A key prediction of our analysis is that 
Dcs/Dq = Tcs/T for completely adamant motors with s = 0. Numerical measurements 
of both the active diffusion constant and the effective temperature have been made by Loi 
et al. [isl for an adamantly motorized semi-flexible polymer melt simulated using molec- 
ular dynamics techniques. Our theoretical prediction is plotted in Fig. [1] along with their 
simulation results. The data points nicely fall on the predicted straight line at low- Tcff 
regime; even at relatively high Tcff (> 4T), only a modest deviation from the slope-1 line 
occurs which slowly increases with rising T^s- The deviation is consistent with our scheme 
being perturbative in the kick size. There are of course also larger error bars in numerical 
experiments for the larger Tgff values. It will be interesting to test our theory for partially 
susceptible motors experimentally. Since D^s is independent of s, for a given motor activity 
A, a higher susceptibility yields a lower T^s and thus a larger slope for the linear relation 
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between Dcs/Dq and Tcs/T. 

Besides this relation our theory suggests T^s should increase quadratically with the kick 
step size for adamant motors, just as is experimentally observed for the dependence on the 
Peclet number found both in sedimentation experiments jsl for active colloidal suspensions 
under gravity, and in a numerical study of an active polymer melt 15 |. 
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FIG. 1: (color online). Rescaled diffusion constant versus rescaled effective temperature for a 
motorized semi-flexible polymer melt. The data points (colored symbols) were converted from 
Fig. 18 of Ref. [l^. The solid line with slope 1 gives the result predicted by the present theory. 



In an early study of the glassy dynamics of an assembly of motorized particles Shen and 



Wolynes [18[ pictured the motors as introducing a modification to the Debye- Waller factors 
of the localized particles. Their non-Hermitian variational approach to finding the steady 
state solution of the master equation turns out to be equivalent to closures for dynamic 
moments leading to an expression of the deviation of the total localization strength a from 
its mechanical value a in terms of the motor properties. On the other hand local mechanical 



feedback within the self-consistent phonon theory [2l| gives back an a for the central particle 
that depends on the a of all its neighbors. Combining these two aspects allows a self- 
consistent determination of mean- field (a, a) solutions and identification of stability behavior 
accordingly. 

Assuming Su = Sd = s, the second moment closure given in reduces to a simple 
expression 



a — a 



a 



* ^ 2 J dD~o 



(8) 



Thus for s = 1/2 {T^s = T) chemical noise does not modify the mechanical stability (a = a); 
for s < 1/2 (Tcff > T) stability is weakened (a < a) whereas for s > 1/2 (Teg < T) stability 
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is enhanced {a > a). 

Here we illustrate the connection to glassy dynamics by carrying out a self-consistent cal- 
culation (see Ref. 



for the detailed procedure) for a motorized version of the cytoskeletal 



network modeled as a "cat's cradle" with excluded volume [16|. Superimposing the phase 
diagrams for active and equilibrium cases over the same parameter ranges as shown in Fig. [2] 
demonstrates that phase boundaries and characteristic densities scale just as we would pre- 
dict from the kick length expansion. 

The phase diagram shows diverse mechanical phases and possible transitions in between 
depending on the density (p) of the constituents and the effective stiffness (P'y) of the 
filamentous elements. In particular the coexistence of a loosely-tethered (LT) mobile phase 
and a repulsive-glass (RG) phase in the multiple-solution (MS) regime allows one to estimate 
the kinetic (p^) and thermodynamic (p^) glass transition densities by matching the free 
energies of the mobile and glassy states. 
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FIG. 2: (color online). Superposed phase diagrams for the equilibrium and active cases for a cat's 
cradle with excluded volume having Dq = 0.025 fim s~ -1, s = 1, K = 10s"\ and / = 0.05 ^m. 
Parameters are comparable to those for a cytoskeleton. The motor kick size / is much smaller than 
the typical crosslink separation which is on the order of l-10/.im. Solid arrows show the changes in 
characteristic behavior, pointing from the phase boundaries and glass transition densities for the 
equilibrium case to those for the active case. Above the upper phase boundary (the crossover line) 
the system is in the tense phase characterized by a smooth crossover from the loosely-tethered 
(LT) phase to the repulsive-glass (RG) phase as the crosslink density p increases. Below the lower 
phase boundary (marked by pth) the mobile phase becomes unstable and only the glassy state can 
be stabilized. Ref. [l^ gives details of the analysis for the equilibrium case. 
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With both susceptible (s = 1) and small-size (/ = 0.05/im) motor kicks, the variety of 
mechanical phases remains intact; but the phase boundaries shift due to the nonequilibrium 
effects. The decrease in leads to an enlarged stability region for the glassy state, while 
increases in glass transition densities and suggest that susceptible motor kicking can 
help resolve local constraints and thus allows a deeper descent into the energy landscape, 
leading to a denser packing where eventually structural rearrangements become too slow to 
be observed. The downward shift of both phase boundaries is simply captured by rescaling 
the equilibrium phase diagram with the inverse temperature ratio (3cs/f3 given by Eq. ([7]). 
Adamant motor kicking gives rise to the opposite effect heating the system above its thermal 
temperature and reducing the glass transition density. Therefore equilibrium phase diagrams 
are still valid in the small kick limit yet with T being replaced by T^s. 

The fact that the master equation precisely reduces to a Fokker-Planck equation allows 
us to apply the effective temperature to analyze the dynamics/kinetics of systems that 
are structurally out-of-equilibrium, i.e., jammed. The slow relaxation of the cytoskeletal 
mechanics 22[ thus is related to glassy behavior in the aging regime. Quantitatively, jammed 
cytoskeletal systems should be described by the aging theory for structural glasses developed 
by Lubchenko and Wolynes [l^]. It is only necessary to replace in their theory the absolute 
ambient temperature which describes vibrations in glasses by the effective temperature of 
the motors. The fictive temperature depends then on the history of the cytoskeleton motor 
activity and preparation. This mapping should make it possible to predict the dependence 
on motor properties of the rheology of active matter. 

Support from the Center for Theoretical Biological Physics sponsored by the NSF (Grant 
PHY-0822283) is gratefully acknowledged. 
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